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The universality of the directed polymer model and the 
analogous KPZ equation is supported by numerical simula- 
tions using non-Gaussian random probability distributions in 
two, three and four dimensions. It is shown that although in 
the non-Gaussian cases the finite size estimates of the energy 
exponents are below the persumed universal values, these es- 
timates increase with the system size, and the further they 
are below the universal values, the higher is their rate of in- 
crease. The results are explained in terms of the efficiency of 
variance reduction during the optimization process. 
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The directed polymer model [1] and the analogous 
KPZ equation [2] have drawn significant attention re- 
cently [3]. The original introduction of the model [1] con- 
sidered directed polymers in random media, and more re- 
cently it has been applied to model other processes such 
as tearing and cracks formation [4]. The study of the 
model is focused on the characteristics of the optimal 
path which connects two distant points in a random en- 
vironment. Passing through any point of that environ- 
ment is associated with a random valued cost, usually 
expressed in units of energy or time. Unlike the simple 
case of constant cost, the path of minimal cost (the opti- 
mal path) is usually not the straight line which connects 
the two points, and its mean transversal distance from 
the straight line which connects its endpoints grows as 
a power law with the distance between the endpoints. 
Another power law behaviour is associated with the re- 
lation between the variance of the energy cost (or time 
cost), and the distance between the endpoints. These 
two power laws are expressed in terms of their expo- 
nents, which are connected by a scaling relation. The 
random environment is usually simulated by a lattice 
whose bonds (or sites) are assigned with random values, 
and it is generally believed that the characteristics of the 
optimal path are universal, i.e. apart from very special 
cases, the exponents depend only on the dimension of the 
lattice and not on the details of the lattice structure, or 
the type of randomness associated with it. 

The universality hypothesis was recently challenged by 
numerical results [5] which showed that for dimensions 
higher than 2, the values of the estimated exponents do 
depend on the type of the probability distribution which 
is associated with the random lattice. In particular, non- 
Gaussian distributions yield lower values for the expo- 



nents compared to Gaussian distribution values, which 
are the persumed universal values. This publication has 
raised strong opposition in the form of three recent ar- 
ticles [6,7,8], which argued that the estimates presented 
in [5] are affected by finite size effects, and thus they are 
not the asymptotic values of the exponents in these cases. 
Two of these articles [6,7] independently suggested that 
the estimates presented in [5] are disturbed by directed 
percolation effects. However, none of these opposing ar- 
ticles presented any alternative numerical estimates to 
those presented in [5] . 

This article presents results of numerical simulations 
for the three and four dimensional cases studied in [5], 
but for larger lattices. The results show that the es- 
timates presented in [5] are indeed below the asymp- 
totic values of the exponents. Moreover, it is shown that 
though convergence to asymptotic values happens at lat- 
tice sizes much bigger than studied, the data suggest the 
presence of the universal values in all cases. In the second 
part of the article the directed percolation explanation 
presented in Refs. [6,7] is challenged, and an alterna- 
tive explanation is presented in terms of the efficiency of 
variance reduction during the optimization process. 

In two dimensions, the directed polymer model can 
be described as follows: In a square lattice which was 
cut along its diagonal and oriented as a triangle whose 
apex is up and the diagonal is its base, we assign to each 
bond or site a random number. Of all the paths leading 
from the apex (the origin) to the base, we refer only to 
those whose direction is always down to the base. For 
each one of these paths we calculate the sum of the ran- 
dom numbers along it, which defines its value. We focus 
our interest on the path of minimal value which we call 
" the optimal path" , and define the space exponent v by 
D t" , where D is the mean distance of the endpoint 
of the optimal path from the center of the base, and t 
is the size of the triangle. In terms of a polymer which 
orients itself along the optimal path, t is the length of the 
directed polymer chain. Another exponent which char- 
acterizes this model is the energy exponent w, defined by 
as ^ i", where E is the the value of the optimal path 
(or the energy of the directed polymer), and ctb is its 
standard deviation. These two exponents are related by 
the scaling relation w = 2i/ — 1 [9] . In the two dimen- 
sional case, the values of the exponents are ~ 2/3 , 
uj ~ 1/3, and the values 2/3 and 1/3 are considered to 
be exact. In dimensions higher then 2, there are only nu- 
merical estimates for the values of the exponents, and in 
the three and four dimensional cases, the accepted values 
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for the exponents are v ~ 0.62, lo ~ 0.24, and v ~ 0.59, 
a; ~ 0.18 respectively [10,11]. However, these numerical 
estimates were obtained for random values taken from 
Gaussian distribution, while non-Gaussian distributions, 
though preserving the space exponents [11], yield lower 
estimates for the energy exponents [5,11]. The depen- 
dence of the estimated energy exponents on the form of 
the probability distribution suggests non-uniersality of 
the model [5], and a break of the scaling relation. It is 
this discrepancy which is adressed in our study. 

Our numerical simulations for d=2 and d=3, were per- 
formed with random bonds whose values were taken from 
uniform distribution. The results are presented in terms 
of local exponents, computed by log4(F(2t)/y(t/2)), 
where V is either D or oe- Since we are interested in 
the relation between the local exponents and the uni- 
versal values, the results are presented in Figure 1 by 
normalized local exponents, where the estimated values 
are divided by the pcrsumcd universal values: v = 2/3, 
u) = 1/3, in 2 dimensions, and v = 0.62, w = 0.24, in 3 
dimensions. The results of Figure 1 can be summarized 
as follows: 

1. All local exponents steadily increase with t. 

2. In both 2 and 3 dimensions the local space expo- 
nents have practically reached their universal val- 
ues. 

3. The local energy exponents in the three dimen- 
sional case are lower than those of the two dimen- 
sional one, but their rate of increase is higher. 

The results related to the two dimensional local energy 

exponents supply a test for the validity of a recent sug- 
gestion proposed in [8]. According to this suggestion, 
the difference between the local exponents and their 
asymptotic values decreases in proportion to t where 
7 = 0.23 ± 0.02. However, the data presented in Figure 

1 suggests that the convergence towards the asymptotic 
value in this case (and in this range of chain lengths) is 
characterizes by logarithmic rather than power law de- 
pendence on the chain length. The probability of power 
law dependence, which is determined by the measure- 
ment errors, was calculated and found to be nil. 

The numerical simulations of the four dimensional 
case were performed on the random site version of the 
model [5], with two types of probablity distributions: 
the uniform distribution, and the arched distribution: 
p{r) ~ (1 — |r|)~^/^, -1 < r < 1. In both cases the 
local space exponents approach their Gaussian values of 
~ 0.59, but in both cases the local energy exponents are 
much lower then their Gaussian values. Moreover, the 
local energy exponents initially decline, and the crucial 
question is whether, as suggested in [5], their minimal 
values are also their asymptotic values. Thus, in Figure 

2 the estimates of the local energy exponents are shifted 



by their minimal values: 0.116 in the uniform distribution 
case, and 0.048 in the arched distribution case, showing 
the difference between the estimates and their minima. 
As can be clearly seen, after the fall comes the rise, and 
the local exponents of the arched distribution case, which 
are lower than the local exponents of the uniform case, 
also increas faster, similar to the picture presented in 
Figure 1. 

In conclusion, Figures 1 and 2 present coherent picture 

of systems which have not yet reached their asymptotic 
local exponents, and whose local energy exponents are 
much more sensitive to the form of the probability distri- 
bution then their local space exponents. The local space 
exponents are close to the Gaussian values, while the fur- 
ther is a local energy exponent below the Gaussian value, 
the higher is its rate of increase. These results provide in 
our opinion a solid support for the universality hypothe- 
sis. 

As mentioned above, both Refs. [6] and [7] suggested 
an explanation of the results of [5] in terms of directed 
percolation effects. In the directed percolation model, the 
random bonds (or sites) are assigned with values taken 
from the bimodal (0,1) distribution, with probability p 
to get zero valued bonds. The set of sites which are 
connected by zero valued bonds to the origin compose a 
directed percolation cluster. In the case of p > Pc, where 
Pc is called the critical probability, there is a finite proba- 
bility to get clusters of infinite length. From the directed 
polymer model point of view, in the case of bimodal (0,1) 
distribution and p > Pc, after an initial search of random 
duration, the optimal path orients itself along a directed 
percolation cluster, and thus the mean energy and energy 
variance of the optimal paths reach constant values which 
are independent of system size, and the energy exponent 
is zero valued. 

In [6] it was shown that in the case of p < Pc, the 
pattern of the local energy exponents is similar to those 
found in the worse cases studied in [5]. Both [6] and [7] 
address the the fact that the higher is the low end of the 
random site continuous distribution, and thus the more 
it resembles the low end of the bimodal distribution, the 
lower are the local energy exponents. Both of them also 
argue that the much lower than universal values recorded 
in [5] for higher dimensions, are the result of the lower pc 
of directed percolation in higher dimensions. A lower Pc 
implies that the same continuous distributions resemble 
more faithfully the critical bimodal distribution with its 
zero valued energy exponent. 

In their reply to the comment of [7] , Newman and Swift 
[12] have already challenged the directed percolation ex- 
planation, and argued that results obtained for the uni- 
form distribution are unlikely to be an outcome of bi- 
modal distribution (or directed percolation) effects. To 
this critique we would like to add another: The crux of 
the matter in the bimodal distribution case, both above 
and below Pc, is the ability of the optimal paths to ori- 
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ent themselves along directed percolation clusters with- 
out any change in their energy values [13]. No such situ- 
ation can occur in the continuous distribution case. This 
fundemental difference is the reason for another differ- 
ence between the bimodal and continuous cases: In the 
bimodal case, not only the energy variance of the opti- 
mal paths grows slowly with lattice size, but so does also 
their mean energy, whose growth rate approaches loga- 
rithmic dependence on lattice size as p approaches Pc- 
Actually, the slow growth rate of the mean energy, com- 
bined with the fact that no optimal path can decrease 
its energy, lead inevitably to the slow growth rate of the 
energy variance in the bimodal case. In contrast, in all 
the continuous cases studied in [5], the mean energy is 
almost proportional to the lattice size. The difference 
between logarithmic dependence and linear dependence 
is a matter of substance, not of a degree. 

Moreover, the whole discussion untill now focused on 
the cases of lower than universal values of the local ex- 
ponents. But opposite cases also exist. Figure 3 presents 
local exponents obtained for a simple two dimensional 
lattice whose bonds are assigned with the negative values 
of random numbers taken from the Gaussian distribution 
and raised to powers of fc = 2, 4, 6. The probability dis- 
tribution in these cases is p{r) ^ \r\^(^^-^)/'^c-^/'2Wf'''-^ 
r < 0, and they all have massive left tails whose length 
grows with k. As can be seen in Figure 3, the longer is 
the left tail, the higher are the local exponents, which are 
even higher then 1/2 at some stages. 

Since the optimal paths orient themselves along lower 
than average energy bonds, the important part of 
the bond probability distribution is its left hand side 
(L.H.S.). The numerical evidence which should be ex- 
plained can be summarized as follows: A decreasing 
(right tailed) L.H.S yields initially lower than universal 
local exponents. An increasing (left tailed) L.H.S. might 
yield initially either lower or higher than universal values, 
the longer the left tail, the higher are the local exponents. 

In the following, an alternative to the directed percola- 
tion explanation is presented. This alternative explana- 
tion relies heavily on the form of the optimal paths energy 
distribution, which for long enough paths arrives eventu- 
ally at a fixed functional form. The persumed universal- 
ity of the model is connected with the assumption that 
apart from very special cases, like the bimodal distribu- 
tion at p > Pc, this functional form is independent of the 
initial probability distribution of the random bonds. The 
asymptotic functional form resembles a Gaussian whose 
left tail is longer than its right tail, and the higher the 
dimension, the more skewed to the left is the form. 

Two different mechanisms determine this form of the 
asymptotic probability distributions: The first is the 
growth of the system, which causes the energy of any 
random path to have eventually a Gaussian distribution, 
regardless of the random bond distribution. The sec- 



ond mechanism is the optimization process, which favors 
lower values, and thus transfers weight from the right side 
to the left side of the optimal paths energy distribution. 

Let us consider a right tailed optimal paths energy dis- 
tribution, and a left tailed one. In both cases a major 
part of the variance comes from the tail. When the opti- 
mization process is performed on the right tailed distri- 
bution, it transfers probabilistic weight from the tail to 
the center, thus reducing effectively the variance of the 
distribution. In contrast, when the optimization process 
is performed on the left tailed distribution, it transfers 
weight from the center to the left tail, thus increasing the 
variance. In the case of a symctric distribution, weight 
is deleted from the right tail and added to the left tail, 
but since the center also moves to the left, the variance 
is reduced. Actiially, even in many cases of pure left 
tailed distributions, the movement of the center to the 
left more then compensates for the additional weight at 
the left tail, and the result is a decrease in the variance. 

As a result of the optimization process, the asymp- 
totic probability distribution is skewed to the left, the 
rate of growth of the variance is lower then 1, and the 
mean energy is not strictly proportional to lattice size. 
Of course, the higher the dimension, the more nearest 
neighbours has each site, and the more efficient is the op- 
timization/variance reduction process. This is the reason 
for both the more skewed to the left asymptotic proba- 
bility distribution, and the lower values of the exponents 
in higher dimensions. 

Moreover, if in the intermediate stage the optimal 
paths distribution is less skewed to the left than its 
asymptotic form, the variance reduction efhciency is 
higher than in the final stage, and thus the local ex- 
ponents are lower than their asymptotic (universal) val- 
ues. The opposite happens in the case of an intermedi- 
ate distribution which is more skewed to the left than 
the asymptotic form. It should only be added that the 
form of the optimal path probability distribution resem- 
bles initially the form of the random bond probability 
distribution, and this analysis explains all the numerical 
evidence summarized above. 

A recent article [14] discusses the random bond proba- 
bility distribution which woTild load to the fastest conver- 
gence of the local energy exponents to their asymptotic 
values. According to the above analysis, the best choise 
is simply a random bond distribution which resembles 
the asymptotic form of the optimal paths distribution. 

The last issue which should be considered is the larger 
deviations from universal values recorded for the same 
random site distributions in higher dimensions. Accord- 
ing to the above analysis, the reason for these larger 
deviations is not the smaller diferrence between these 
distributions and the critical bimodal distribution, but 
the bigger difference between these distributions and the 
more skewed to the left asymptotic form of the optimal 
paths probability distribution. 
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In conclusion, this study supports the universality 
of the directed polymer model in two ways. First, it 
is shown that finite size effects do influence the esti- 
mated local exponents of the non-Gaussian distributions 
in higher dimensions. Then, these finite size effects are 
explained by the form of the optimal paths probability 
distribution, and its influence on the efficiency of variance 
reduction during the optimization process. 

We wish to thank the German Israel Foundation (GIF) 
for financial support. 
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FIGURES 

Fig. 1. The local exponents of v and lo (normalized 
to the universal values), as a function of the chain length 
for d=2 and d=3. The random bonds values are taken 
from uniform distribution. 

Fig. 2. The local exponents of a; as a function of the 
chain length for the random site lattice in d=:4. The re- 
sults are shifted by their minimal values: 0.116 for the 
uniform distribution, and 0.048 for the arched distribu- 
tion, thus presenting the the difference between the local 
exponents and their minima. 

Fig. 3. The local exponents of a; as a function of the 
chain length. The results are shown for 3 random bond 
probability distributions in d=2. The bonds are assigned 
with the negative values of random numbers taken from 
the Gaussian distribution and raised to powers oi k = 
2,4,6. 
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